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Nucleation and propagation of a shear fault is known to be the result of interaction and coalescence of many 
microcracks. Yet the character and rate of the microcracks' interactions, and their dependence on the 
three-dimensional stress state are poorly understood. Here we investigate formation of microcracks during 
sandstone faulting under 3D-polyaxial stress fields by analyzing multi-stationary acoustic waveforms. We 
show that in a true three-dimensional stress state (a) faulting forms in a orthorhombic pattern, and (b) the 
emitted acoustic waveforms from microcracking carry a shorter rapid slip phase. The later is associated with 
microcracking that dominantly develops parallel to the minimum stress direction. Our results imply that 
due to inducing the micro-anticracks, the three-dimensional (3D) stress state can quicken dynamic 
weakening and rupture propagation by a factor of two relatively to simpler stress states. The results suggest a 
new nucleation mechanism of 3D-faulting with implications for earthquakes' instabilities, as well as the 
understanding of avalanches associated with dislocations. 

The mechanisms of rock faulting were mainly addressed by two macro-mechanical models, namely ander- 
sonian 1 and slip models 2 that are compatible with field and experimental observations 3 5 . While the classic 
andersonian model (based on the Coulomb-Mohr criterion) predicts bimodal faulting (conjugate set of two 
faults), the slip model explains the polymodal fault patterns (orthorhombic sets of four faults) 2,6,7 , observed under 
a true triaxial stress state (ai > rj 2 > ff 3 ). From a micro-mechanical point of view, based on acoustic emission 
patterns in conventional triaxial tests (i.e., a 1 > a 2 = ff 3 ), the shear faulting is related to the interaction and 
coalescence of brittle, dilational microcracks 810 . The prevailing observation is that these microcracks are pref- 
erentially aligned parallel to a l — a 2 plane, namely, normal to the least compressive stress 6,8,10 . On the other hand, 
the microcracks that are associated with deep-focus earthquakes, which occur under high pressure and high 
temperatures, are possibly oriented perpendicular to maximum compressive stress, and termed anticracks 11,12 . An 
important feature of deep-earthquakes, which may be related to their anticrack source, is that they start up 
significantly faster (i.e., shorter rise time) than do intermediate or shallow depth earthquakes 13,14 . Here, we show 
that shorter rising time is not restricted to deep-focus earthquakes, and is observed in unstable slip under a true- 
triaxial stress state. We monitored the detailed time evolution of 3D faulting that led to polymodal, complex fault 
systems, and visualize this dynamic evolution of rupture-forming faults under true-triaxial test conditions (TTT). 
We found that the nature of unstable slip under TTT can be different from unstable slip under conventional 
triaxial tests (CTTs). We show that fault-nucleation in TTTs is associated with shorter generic phases of emitted 
waveforms which imprints unique signatures in appropriate phase spaces, and they differ from rupture fronts in 
CTTs. Our experiments further highlight a potential new weakening mechanism which implies that micro faults 
reach about a factor of two faster to steady state frictional strength. 

Results 

Our experimental system is schematically shown in Fig. la. A cubic sample of Fontainebleau sandstone (80 X 80 
X 80 mm) with 3 acoustic emission (AE) transducers per face is loaded normal to its faces. The stresses a 2 = 
35 MPa and a 3 = 5 MPa are maintained constant during the experiment (see section 1 and 4 of supplementary 
materials for other test types). The superimposed best-located AE events on the reconstructed micro X-ray 
patterns showed that the evolution of faulting follows "diffusive" fracture patterns, forming the final shape of 
fault patches (Figs, lc, d, SI). This is a typical example of composite self-organized cracks, observed frequently in 
cylindrical CTTs 15,16 . Despite this similarity, we found that the evolution of acoustic signals corresponding to 
rupture fronts is different in TTTs. 

To explore the above observations, we used network theory modeling on the recorded multi-stationary acoustic 
signals ([17]). We use three parameters of this theory to characterize the AE distribution and evolution. The Q 
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Figure 1 | Study of Orthorhombic faulting under the true triaxial tests, (a) a schematic of the experiment procedure and remote stress field configuration 
have been shown where the numbers on each face of the cube shows the piezoelectric transducers, (b) The evolution of the main driving stress field (oj) 
and accumulated recorded acoustic emissions during a few hours of the experiment. The stairlike trend of the loading curve is due to pausing the test in 
order to measure permeability in three directions. (c,d) Superimposed events on the best-located events superimposed on the reconstructed fault patterns 
from X-ray images and the temporal/spatial evolution of the events (i.e., event-sets' sequences). The first two event-sets show the slow growth of the 
external fault-system's patches while the fast-growth of the internal branches of the Orthorhombic faulting system is observed around the main stress 
drop (Fig.Sl). The numbers shows 4 main fault sets, (d) A close view of 3D - fault system. 



parameter is the 'modularity' of the network that characterizes the 
number of modules of network of the AE event. The B.C parameter is 
defined as the "betweeness centrality" that approximately is corre- 
lated with the magnitude of the AE events. The lambda parameter is 
defined as the maximum eigenvalue of Laplacian of the network (see 
Methods part for definition and details), which scales with the 
inverse of Q parameter in our experiments. Applying these para- 
meters on acoustic signals, we found three generic time scales, cor- 
responding to the following phases in Q-profiles or modularity of 
corresponding functional networks 17 : (1) a nucleation and main 
deformation phase in the form of initial strengthening (2) a fast-slip 
or weakening phase, and (3) a slow slip or decelerating stage (Fig. 2a). 
Comparing the evolutionary phases of micro-cracks to dynamic 
stress change from macro-sliding events (recorded by strain gauges) 
reveals similar time characteristics and trends (Fig. 2c, d). In using 



the reciprocal of Q-profiles, we magnify the initial strengthening 
phase (we call it R-profile) to obtain the maximum strength (Fig. 2b). 

We found (SI, Figs. S2) that the duration of the fast-slipping phase 
of precursor events does not effectively change under a wide range of 
loading conditions including different strain rates, acoustic signal 
feed-back control test, smooth and rough fault frictional interface 
with 150 and 200 MPa confinement pressures (Westerly granite 
samples), and double-direct shear test with quartz-gouge materials. 
These observations strongly indicate that the fast-slip phase is a 
generic phase in micro-cracks, which is not sensitive to variations 
of loading conditions in CTT or conventional friction tests. We note 
that a nearly invariant time scale of the second phase has been 
recently reported 18 on engineering polymeric materials (and inter- 
faces) and has been shown to be insensitive to the dimension of 
loaded samples. Comparing the fast-slip times of micro-cracks under 
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Figure 2 | Q-profiles resemble dynamic stress changes in micro-scale failures, (a) The main phases of Q-profiles includes three main stages: (1) initial 
strengthening, (2)fast-slip or weakening phase and (3) deescalating or slow-slip phase. The profile is from CTT- westerly granite experiment 16 , (b) The 
maximum dynamic strength is inferred from R-profile as the inverse of Q-profile which places emphasis on the evolution of the first phase. We have 
shown —20 acoustic excitations from TTT (FTB4). The inset illustrates a simplified R-profile. The weakening stage can be described by the time- 
weakening model which prescribes that strength of heterogeneity weakens during a characteristic time (tu). (c, d) Dynamic stress changes from 5-macro 
stick-slip (Saw-cut fault evolution 16,17 in Westerly Granite) resemble similar evolutionary phases and time-characteristics with Q-profiles in micro- 
ruptures (d). 



CTT, from different rocks (with dominant silicate minerals), con- 
crete samples(Portland cement), thermally treated Westerly granite 
specimens (850°C and 250°C), Polymers and steel alloy implied that 
this generic phase is linked to the ductility of materials, and changes 
in energy dissipation mechanisms could alter the duration of the 
second phase (section3 of SI and Fig. S.3). For instance, this period 
is longer for concrete -as cement based materials-due to the cement 
component and Calcium-Silicate-Hydrate (C-S-H) structures 
(Fig.S.3c). This is due to additional sources of energy dissipation 
rather than the creation of new surfaces 19 . For the same reason ther- 
mally treated (850°C) westerly granite specimens, characterized by 
much higher intergranular and intragranular crack density than 
250 U C specimens, -shows a fast-slip time of —40 (is in comparison 
to the later specimen (—20 us). Also we discover that this phase is 
about —70 us for Poly (tetrafluoroethene), close to —60 us for 
PMMA 18 . This finding suggests that longer fast-slip phase likely 
indicates a signature of plasticity in micro-cracks (i.e., significant 
ductile cracking). 

Interestingly, examining the same interval for a true triaxial test 
shows a shorter time interval (—10 us) for the same sandstone 
(Fig. 3a (and inset), b). We have observed the mentioned phenom- 
enon from different events, occurring in different times and positions. 



This observation implies a different pulse shape for acoustic- cracking 
noises under 3D boundary conditions while a universal left-asym- 
metric shape of the pulse is maintained for both CTT and TTT 
(section3 of SI). To further analyze the shorter fast-slip phase, we 
use a simple effective temperature model ([18, 20] -section2 of sup- 
plementary materials). The model estimates that the cracks under the 
triaxial stress field carry a smaller amount of energy (at least 40% of 
CTT), indicating different nature of cracking. We support the idea of 
micro anti-cracks with observation of significant activity of micro- 
cracks in the direction of the minimum stress state, forming anit- 
cracks (Fig. 3c-sections 4,5 of SI - Fig.S.l,S.6, S.18-19). We speculate 
that anti-cracks form rupture fronts, induced by far field stress state, 
results in a faster slipping or weakening period. 

To see whether micro-anticrackings changes the shape of wave- 
forms, we examined typical waveforms during the formation of 
anticracks in a multi-anvil test (High pressure -High temperature: 
HP-HT) on the Olivine samples (section 5 of SI - Fig.S23-S24), 
which revealed —8 us of the fast-slip phase. Thus, this observation 
of events from the HP-HT experiment confirms the shorter rising 
time of natural deep-focus earthquakes 1314 ; a feature that is similar to 
the shorter weakening phase of rupture fronts in TTTs. Interactions 
of such anticracks with each other likely form the final fault surface 
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Figure 3 | Q-profiles and fast-slip phase in TTT and CTT. (a) Q-profiles from a TTT test (FTB4) show an —8-10 us duration of the fast-slipping phase 
(inset: 150 events from TTT-FTB2; see section 4 in the supplementary information). (b)The same sandstone under the conventional test shows a longer 
fast slip phase (—23-25 us). Inset shows 45 events from the CTT. Q(f 0 ) is the initial or rest value of modularity, (c) Visible (anti-) cracks parallel to 
minimum stress field in TTT. 



which can be explained with a microcrack interaction model such as 
the three-dimensional solution of Eshelby 6 . The concept of anti- 
cracks previously has been also used to analyze the failure process 
in high porosity sandstone as a theory of compaction bands 21 , to 
characterize the triggering mechanism of snow avalanches 22 as well 
as the localized dissolution of limestone 23 . We conclude that anti- 
crack form of rupture fronts and possible substantially lower fracture 
surface energy are the main mechanisms of our observations. Due to 
the frequent observation of micro-faults with dominant double cou- 
ple source mechanism (i.e., shear component), we propose that true 
3D stress states can induce rapid weakening on 3D-faults, faster than 
conventional weakening mechanism(s). This fast-release of energy is 
associated with fast-acceleration 24 which shortens the slip-weaken- 
ing distance. The nature of this slip-acceleration is different from one 
demonstrated in 24 which is associated with fault wearing and gouge 
formation; here we assign the quickened phase to grow of micro- 
anticracks. 

Furthermore, we show that network parameter spaces of TTTs' 
events, collectively, are different from CTTs' events. Considering 
mean-temporal of maximum Laplacian (l max ) an d spatio-temporal 
mean of centrality (log <B.C>) on networks (see Methods part for 



definition), we find that recorded events from TTTs and CTTs col- 
lapse in /i max — log <B.C> 1 plane (Fig. 4a), which confirms recent 
similar observations from the rock-frictional interfaces ([17]- 

Fig.S8,Fig.S.14 and Fig.S.17). Considering that Q<xA max 1 and 
Qccer -1 (cr is the temporal average of stress change during the mon- 
itoring interval), we deduce that events with higher 2 max indicate 
ruptures with higher weakening rate and then high rupture velocity 34 . 
As we have shown in Fig. 4a-inset, the conventional tests imprint 

three main clusters of events in A max — log <B.C> , correspond- 
ing to the energy spectrum of ruptures 17 . Comparing 

/i max — log <B.C> 1 plane on TTT and CTT (Fig. 4c), it turns 
out that in TTT, we cannot recognize distinctive clusters, indicating 
that probability of finding "abnormal" events in terms of the rup- 
ture's speed regimes under true triaxial conditions is rare. Thus, these 
observations imply that energy of events is more localized in TTTs. 
This indicates that the shapes of waveforms are also bounded rather 
than CTT cracking noises (Fig.S.15-18). 

To complete our analysis, we speculate that the widely discussed 
failure criterion of materials under TTTs 2,25-27 can be inspected in 
the scale of tiny- amplified events (i.e., micro-failure criterion). To 
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Figure 4 | Global and local networks' parameter spaces from precursors failures on TTT and CTT. (a) A max — log <B.C> 1 parameter space on —4500 
events (including weak events) an hour prior to the main-final faulting development in FTB4 (section 3 of SI). Inset shows events from the CTT at the 
same plane with three main classified classes corresponding to three classes of rupture fronts energies 17 . Class I, II and III corresponds with slow, regular 
and very fast weakening rates, respectively. Generally A max oc Q and Qoc c~ 1 (b) a study of independent layers of networks. The sub-network's 
parameter space: collapsing —520 strong events from the FTB4 test in the ] t — J 2 parameter space, where J 1 =R[+R 2 + R ! 3 and ] 2 = R[ R ! 2 + R 2 Rl + R[ R 3 ■ A 
power law relation satisfies an excellent collapsing of events in the Jj — J 2 phase diagram as: J 2 oc /j\&~2.6. This scaling confirms Reches' proposed failure 

criterion 2 . In the inset, we have shown the R[ — Yr plane with Yr = y (R[ —R' 2 ) 2 + (R^—R^) 2 + (R 3 —R[) 2 -proportional to z act - for the same events 
(Also see Fig.S.5-Fig.S.15). (c) a collection of —3000 events from CTTs and TTTs (see the list of CTTs in SI) in the normalized Q —log <B.C> 
map. 



envisage such micro-failure criterion, we consider the possible rela- 
tion between independent layers of functional multiplex acoustic 
networks (Methods part), in which an independent network of sen- 
sors is assigned for each main direction(X,Y and Z). Then, for each 
occurred rupture front and for each sub-network (layer), we calculate 
the reciprocal of the minimum modularity index, namely "maximum 
resistivity" indicates the maximum resistance against motion 
from i th layer or sub-network in the first phase of Q-profries - 
Fig.S7 and Fig.S9-S17). Since R\ occurs in the main deformation 
phase of Q-profiles (i.e., first stage), then it is proportional with the 
maximum strength in that direction (Kjao J ^' r ' z '). Remarkably, 
the events forming the Polymodal fault-system collapse in the / t — J 2 
parameter space in which h = R[ +R' 2 +R' 3 and J 2 = R[R 2 + R 2 R' 3 + 
R'iR^ (Fig. 4b), represent an excellent power-law: } 2 oc ]\,b~2.(s (also 
see section3 in SI). Indeed, Jx and } 2 are proportional with the stress 
invariants where J l oc a l + a 2 + o" 3 and ] 2 oc a 1 a 2 + a 2 a 3 + <7ia 3 . 
Interestingly, the obtained power-law and similar relations (Fig. 3b- 
inset) on other parameters match with the macro-failure criterion 



presented by Reches 2 , Mogi 25 and Haimson 26 , indicating that possibly 
a similar mechanism governs both precursor and final failures (i.e., 
significant drop of the driving stress). 

In addition, based on Mogi's approach 25 , under TTTs, the ductility 
of rocks decreases. To show the validation of this character in micro- 
failures, we note that the signature of ductility (plasticity) is encoded 
in weakening time (SI, Fig.S3), so that it generally increases the 
duration of this phase which leads to the further decrease of l max 
and the possible formation of distinct cluster. Then, due to the 
absence of the lower cluster (cluster I in Fig. 4a-inset; Fig. 4c) in 

A max — log <B.C> from TTTs, we speculate that a similar 
approach holds in micro-failures. The same logic can be followed 
for a i — z oct plane which is proportional with R[ — Yr (Fig. 4b-inset) 

in which i ort . = - ((o^ - a 2 f + (<r 2 - 17 3 ) 2 + {a x -ff 3 ) 2 ) ' 2 - In other 

words, we have shown that the long-standing problem of failure 
criterion in a 3D stress field has its origin in each single recorded 
precursor event (i.e., local failure). 
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Discussion 

Our results have significant implications for earthquake seismology, 
rock-mass stability and fluid migration and 3D permeability aniso- 
tropy in fractured rocks. Our results provide a link to prove that 
accommodation of 3D strain due to polymodal fracture patterns 
change strength (encoded in failure criterion). For the first time, 
we showed that micro-anticracking can change slip-weakening rate 
and accelerate the rate of energy release. This implies a sharper 
source time function for further modeling of our experiments. 
Indeed, our results can be assumed as an extension to detachment 
fronts in micro-faults 18 to a general concept of anti rupture fronts. In 
addition, we showed that deep earthquakes and fault instabilities 
under 3D-polyaxial stress fields diverge little in terms of the stress 
drop invariance (i.e., fast-slip duration), promoting the recent 
approach on similarity of deep earthquakes to their shallow counter- 
parts 33 . However, this similarity does not imply that the same physics 
governs both weakening courses. In deep-focus earthquakes, it is 
likely the nucleation is associated with dislocation density and com- 
plex motion of dislocations is encoded in fast-slip regime (Fig.S.24) 35 . 
Regarding nucleation of anti-cracks in our TTTs, since the grains of 
sandstone can arrange in a tightly packed configuration, the volu- 
metric collapse and cohesion along the crack faces significantly 
shorten the critical displacement 22 . Our results confirm that the 
growth of anti-cracks (as the source of true compaction bands) is 
due to rapid propagation as suggested in 21 . From another perspective, 
studies in atomic friction anisotropy revealed differences in acoustic 
excitations can be related to different crystallographic directions 
where in a certain direction the energy dissipation is minimum 36 . 
Our findings using networking of anti micro-cracking's acoustic 
excitations in TTTs coincides with such observations. The presented 
results might be applicable for emitted events from dislocations in 
atomic scales which radiate elastic energy in terms of acoustic emis- 
sions 28 . In terms of the analysis of micro damages, our algorithms 
provide unique tools to understand the details and physics of multi- 
stationary acoustic waveforms. In fact, these algorithms can be 
employed to extract further information on other sources of acoustic 
excitations. Developing models based on multiplex-functional net- 
works and extending the proposed network analysis to natural earth- 
quakes will be a future project. 

Methods 

Summary of experimental procedures. Cubic-samples of the saturated 
Fontainebleau sandstone are loaded at the University of Toronto, in a custom-made 
true triaxial cell 29 . The triaxial cell can be used to measure permeability independently 
in different directions under raised temperatures. Special attention was paid to keep 
the balance of loads in parallel directions to avoid possible torque. The strain rate for 
all axial directions was 3 X 10" 6 s~\ The waveforms were recorded using 18 
piezoelectric transducers at a 10 MHz sampling rate while triggered sensors are 
employed to detect the events. Another system simultaneously and continuously 
recorded all events from all transducers during the aforementioned experiments. For 
conventional compressive tests, we used the results of previously reported 
experiments 30 while the samples were saturated and the events were chosen from 
precursor rupture fronts. A Fontainebleau sandstone specimen {length — 88 mm, 
diameter — 40 mm) was deformed inside a triaxial. A network of piezoceramic 
transducers (PZT) was used inlO MHz sampling rate with similar criterion to the 
Polyaxial test to detect the events. In the supplementary document (sectionl and 5), 
we have -also-included the summary of other tests (and their multiplex networks 
attributes) on Westerly Granite samples under different loading conditions and 
geometry as well as smooth and rough rock-frictional interfaces' 6 ' 17 . Moreover, we 
used our ultrasonic imaging triaxial cell to record the events emitted during testing of 
850°C and 250°C thermally treated Westerly granite specimens (125 mm long and 
50 mm in diameter) tested under 25 MPa of confining pressure. In addition, we 
analysed multi-stationary waveforms from a Multi-anvil test (High pressure -High 
temperature tests) on olivine samples (section 5-SI). 

Functional multiplex networks on acoustic emission waveforms. We define a 
multiplex network on the cubic geometry over acoustic transducers where it includes 
three main layers {the layers correspond to X,Y and Z directions as the minimum, 
intermediate and maximum principle stress). To build networks, we use the 
algorithm introduced in 17 . We built a network on each time step in which the nodes 
correspond with acoustic sensors where elastic "excitement" induces voltage- 
fluctuations in them. Please note that the "functional" friction networks are not 



structural networks; however underlying functionality of constructed networks 
correlates with main mechanical features of crack instabilities (see SI for more 
information and details of algorithms). Two cases were analyzed through this 
research: global networks or interdependent layers and local-sub networks as 
independent layers of mathematical graphs. Each node in functional networks was 
characterized by its degree For a given network with N nodes, the degree of the 
node and Laplacian of the connectivity matrix were defined by 

N 

kj — a^; Ly — ciy — kjSjj where k^a^Ly are the degree of i th node, elements of a 
symmetric adjacency matrix, and the network Laplacian matrix, respectively. The 

N 

eigenvalues are given by^i^"' =A„ ( j4f ) ,m which istheith eigenvector of 

7=1 

the Laplacian matrix (x — 1, TV) .We define A max — { — A N ) as the maximum 
eigenvalue of the Laplacian of the network. The betweenness centrality (B.C) of a 
node is defined as 31 : 

1 N PjS 

BQ= (N-1KN-2) 5 t (1) 

in which p^ is the number of the shortest path between h and j, and p^j is the number 
of the shortest path between h and j that passes i. The spatio-temporal average of B.C 
is indicated by log < B. C > where < . . . > and bar-sign correspond to the spatial (i.e., 
nodes) and temporal averages, respectively. In 17 , we speculate that events from saw 
and rough fault cut -experiments flow in A ma x — log <73.C> 1 plane with three 
distinguished trends, corresponding to rupture regime. 

The Q value per each network {i.e., in each time-step) as the networks' modularity 
characteristic is given by 31 ' 32 : 




in which N M is the number of modules (clusters), L — - kj, l s is the number of links 
in module and d s — k\ {the sum of nodes degrees in module s). Using an 

optimization algorithm (Louvain algorithm- [31]), the cluster with maximum mod- 
ularity (Q) is detected. The temporal evolution of Q values in a few hundred micro- 
seconds evolution of waveforms provides Q-profiles. 

As we have shown in Fig. 2. and Fig.S4, the main time characteristics of the 
evolutionary phases of Q-profiles and their trends from micro cracks coincides with 
dynamic stress profiles in macro -ruptures (a few centimeter scale). Then, we infer 
that the maximum Q (t) corresponds with the maximum stress drop (^ min _) in sub- 
micron scales. The initial drop of Q-profile (as the nucleation and strengthening 
phase) before the weakening phase corresponds with maximum strength of hetero- 
geneity (or asperity). The minimum value of Q(t) (or maximum value of R(t)) 
coincides with *r max . Inferring this 0" max . in each main directions (such as X-Y-Z) 
provides the ingredients to infer local failure criterion. For this reason, we address the 
maximum values of R{t) in three main directions in our TTTs. We also note that for 

all our data set, we have: Qoc Regarding interpretation of log <B.C> \ we 

note that minimum value of log<B.C> generally coincides with the large amplitudes 
of waveforms (Fig.S.6). Then, we infer rupture with high energy will have larger 

log <B.C> 1 . There are some indications that critical temporal curvatures of Q (t) 
represent significant fluctuation in log <B.C> 1 (Fig.S.5). With using the layers of 
multiplex networks and distribution of log <B.C> 1 in each layer (or direction), we 
can infer the tendency of rupture per each recorded micro-event (Fig.S.7). 
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